Astronomy <fe Astrophysics manuscript no. h4352'3 


February 2, 2008 


(DOI: will be inserted by hand later) 





o 
o 



W Hya : molecular inventory by ISO-SWS * 

K. Justtanont 1 , T. de Jong 2,3 , A.G.G.M. Tielens 4 ' 5 , H. Feuchtgruber 6 , and L.B.F.M. Waters 3 

1 Stockholm Observatory, AlbaNova, Dept. of Astronomy, SE-106 91 Stockholm, Sweden 

2 SRON-Utrecht, Sorbonnelaan 2, NL-3584 CA Utrecht, The Netherlands 
<^- ~< 3 University of Amsterdam, Kruislaan 403, NL-1098 SJ Amsterdam, The Netherlands 
Q ■ 4 Kapteyn Institute, P.O. Box 800, NL-9700 AV Groningen, The Netherlands 
Q [ 5 SRON-Groningen, P.O. Box 800, NL-9700 AV Groningen, The Netherlands 

£N| ■ 6 Max-Planck Institut fuer extraterrestrische Physik, Giessenbachstrasse, D-85740 Garching, Germany 

X> : 

, Received ; accepted 

Abstract. Infrared spectroscopy is a powerful tool to probe the inventory of solid state and molecular species in 
I ■ circumstellar ejecta. Here we analyse the infrared spectrum of the Asymptotic Giant Branch star W Hya, obtained 

by the Short and Long Wavelength Spectrometers on board of the Infrared Satellite Observatory. These spectra 
, show evidence for the presence of amorphous silicates, aluminum oxide, and magnesium-iron oxide grains. We 

^ ' have modelled the spectral energy distribution using laboratory measured optical properties of these compounds 

and derive a total dust mass loss rate of 3 x 10~ 10 Mq yr _1 . We find no satisfactory fit to the I3(im dust emission 
feature and the identification of its carrier is still an open issue. We have also modelled the molecular absorption 
bands due to H2O, OH, CO, CO2, SiO, and SO2 and estimated the excitation temperatures for different bands 
which range from 300 to 3 000K. It is clear that different molecules giving rise to these absorption bands originate 
from different gas layers. We present and analyse high resolution Fabry-Perot spectra of the three CO2 bands in 
the 15/xm region. In these data, the bands are resolved into individual Q-lines in emission, which allows the direct 
determination of the excitation temperature and column density of the emitting gas. This reveals the presence of 
«4 ' a warm (~ 450K) extended layer of CO2, somewhere between the photosphere and the dust formation zone. The 

gas in this layer is cooler than the 1 000K CO2 gas responsible for the low-resolution absorption bands at 4.25 
and 15/im. The rotational and vibrational excitation temperatures derived from the individual Q-branch lines of 
'_h ' CO2 are different (~ 450K and 150K, respectively) so that the CO2 level population is not in LTE. 

C/3 '■ 
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H , 1. Introduction weak CO millimeter observations (Wannier & Sahai 

■ - - ■ Bujarrabal et al. 119891 Cernicharo et al. [1997). It is clas- 

One of the characteristics of Asymptotic Giant Branch gified ag & semiregular variable {SRa) from its light curve . 

(AGB) stars is their mass loss. As a star evolves up the The vigual magnitude ranges from 6 to 10 over a pe _ 

AGB, its mass loss rate increases. Just before it evolves off riod of 382 dayg (Lebzelter et aL EJjjjjj. It exhibits Si0 

to become a planetary nebula, it sheds most of its mass (e g ^ Schwartz et aL Qggg clark et aL H2 q ( e .g., 

in a short time, i.e., superwind. Signatures of mass loss Spencer et ^(jgjg Lane et aL [J^ and qh main line 

can be studied in the infrared by taring the dust (e.g., masers (e g ^ Szymczak et jj.^ggg Chapman et al. 11553 . 

BedijnUMTl Justtanont & Tielens |T992) , and m the sub- ft has ft wgak silkate dugt emission at l0 ^ m along with a 

millimetre by using mokcular rotational Imessuch as CO smaUer emission feature at 13fim xhe carrier of the latter 

(ag -' ^ P & MOTris ™ GroenewegenHnSI Skinner ^ gtm bdng debated w Hya ig Qne of ^ firgt AGB starg 

et al. | 1999 ||. £ or w ] 1 j c j 1 JJ2O rotational emission was reported after the 

An example of an early AGB star is W Hya. It la unch o f the Infrared Sp ace Ob servatory (ISO, Kes sler et 

has a relatively low mass loss rate as deduced from the a l . IT9"9"r)|l . Neufeld et al. ifTTMl) and Barlow et al. fTM^ 

reported that these lines were detected in Short and Long 

Send offprint requests to: K. Justtanont w n c , , /owe* j t^tc<\ j.- 1 

* _ , . . Tr ,^ „ . , Wavelength Spectrometers (SWS and LWS), respectively. 

Based on observations with ISO, an ESA project with m- 



. . , , , , . ,, . ,,, Both authors presented models accounting for the emis- 

strumcnts funded by LSA Member States (especially the PI . . 

, . „ r , AT , , , TT ., , sion. Recently, Zubko & Zuckerman (2000) also reported 

countries: trance, Germany, the Netherlands and the United ■" ' » ^ 

Kingdom) with the participation of ISAS and NASA. a model predicting H 2 line fluxes observed by ISO. Also, 
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W Hya has been detected in the ortho ground state line by 
SWAS at 557GHz (Harwit & Bergin I2UU2)> In our paper, 
we mainly discuss the SWS spectra taken during the ISO 
mission. We also present the SWS Fabry-Perot spectrum 
of W Hya centered on C0 2 bands. This band has been 
reported by Ryde I1999|l and Justtanont et al. (p.998). 

2. Observations 

Spectroscopic observations of W Hya were taken aboard 
ISO using both the SWS (de Graauw et al. ITMfJjl and 
LWS (Clegg et al fiMty . The SWS data were taken using 
the standard AOT1 (full scan from 2.38-45.2/im) speed 

3, giving an average resolution of 510. We also obtained 
the spectrum using the full grating resolution of 2000 be- 
tween 12-16/im (AOT6) in order to study the CO2 bands. 
The data were reduced using the data reduction package 
provided by the SWS calibration team based at SRON- 
Groningen. We were able to observe the three CO2 bands 
using the SWS Fabry-Perot (FP) centering on the three 
Q-branches in order to resolve individual transitions. Care 
was taken to deglitch the spectrum due to cosmic ray hits 
on the detectors at the time of observation. As the data 
were taken in small wavelength intervals, each subband 
spectrum had to be stitched together. The final spectra 
were rebinned to a resolution of 410 , with an oversam- 
pling of 6. Also the flux level for each wavelength band 
has been normalised to the AOT6 flux. 

The LWS data were taken from the archive (PI: M. 
Barlow). The reduction was done using the ISAP data 
reduction package provided by IPAC. 

3. Dust and gas mass loss rates 

In this section, we present a model fit to the infrared spec- 
tral energy distribution (SED) in order to derive the dust 
mass loss rate. This is then used to calculate the rate of 
gas which is driven by the dust to the observed termi- 
nal velocity by momentum transfer between gas and dust. 
At the same time, the temperature structure of the cir- 
cumstellar envelope is calculated. The main gas heating is 
by the dust drag while cooling is done through adiabatic 
expansion of the gas and radiative cooling via rotational 
molecular line emission. We also include vibrational cool- 
ing by H2 (Goldreich & Scoville 119761 Justtanont et al. 
I1994J) . Hence the dust and gas mass loss rates are derived 
in a self-consistent manner. The issue of mass loss rates 
is readdressed here as there is a wide range of reported 
rates from 3 10~ 5 (Neufeld et al. IT§5Bjl to 5.2 10~ 8 M Q 
yr -1 (Wannier & Sahai 1986); see also Zubko & Elitzur 
( 2000 ) for a recent reassessment of this controversy. 

3.1. The dusty envelope 

To derive the dust mass loss rate, we have fitted the SED 
of W Hya using ISO SWS and LWS data, together with 
photometry from ground based observations (Wilson et 
al. I1972[) and the IRAS point source fluxes. We adjusted 



the flux levels for the LWS spectrum to coincide with the 
IRAS 60/im flux which also joins the SWS spectrum very 
well. There are broad features present in the LWS spec- 
trum which are due to the internal reflection within the 
instrument at 54 and llO^m. Here, we calculate the radia- 
tive transfer through the dust shell, using a code based on 
Haisch (1979) which takes into account thermal emission 
and reemission and multiple scattering. It also allows for 
a grain size distribution which we assumed to be the same 
as the interstellar grain sizes (Mathis et al. 119771 hereafter 
MRN). We also employ multiple grain species. 

The input parameters are listed in Table ^ The dis- 
tance is taken from the Hipparcos catalogue (ESA lH)9"?|l . 
From the study by Cami l)2()02|l . the emission features be- 
tween 10-20/im are a combination of amorphous silicates 
giving rise to the 10/mi emission, compact AI2O3 which 
gives a feature at 11/zm, spinel (MgA^O^ is responsible 
for the 13/zm feature (Posch et al. 1992J), and the 19/im 
emission which is attributed to Mgo.1Feo.9O. In our at- 
tempt to fit these dust features, we used silicate dust prop- 
erties employed by Justtanont & Tielens l|1992|l . For other 
compounds, we extracted the optical constants from the 
University of Jena database and calculate the extinction 
coefficients using the Mie theory. 

The fit to the SED for a three-dust component model 
is reasonably good for the complete SWS and LWS data 
coverage (Fig. The fit is very satisfactory for all ob- 
served emission features, providing support for the pres- 
ence of amorphous silicates, AI2O3 and MgFeO, except for 
the 13^m feature. A model with spinel which produces a 
fit to the 13/Ltm emission feature predicts an additional 
feature at 16/im, not seen in the observation (Fig. EJ. 
Hence, while spinel provides a good fit to the 13/im band, 
this identification is suspect. Further laboratory studies 
on related materials are required to test this suggestion. 
Onaka et al. ( 1989) suggested aluminum oxide as a carrier 
for the 13/j.m dust emission and following that, Kozasa & 
Sogawa |1908J) were able to explain the 13/um feature us- 
ing crystalline aluminum oxide core-silicate mantle grains. 
However, the optical constants used for aluminum oxide 
were coarse (l/um spacing, Toon et al. I1976JI and the in- 
terpolation of them close to the resonance could result in 
a misleading profile. A possible carrier proposed by Speck 
(200tJ) is Si02 which has features at 9.5, 13 and 20/im 
while Begemann et al. (^9?) suggested a silicate-based 
material. For the moment, the issue of the dust species re- 
sponsible for the observed 13/mi feature is still unresolved. 

At a distance of 115pc, we derive a total dust mass loss 
rate of 3 10 -10 M yr _1 for the three species of silicates, 
AI2O3 and MgFeO, assuming a constant outflow velocity 
of 6 km s _1 (Bujarrabal et al. I1989|l . We note here that 
we assume that the dust and gas velocity is the same as 
the fitting of the SED allows us to derive M^/wdust- We 
will discuss the effect of the dust drift velocity in Sect. 3.2. 
Our input values are very similar to those used by Zubko 
& Elitzur ( 2000 ) but we derive a dust mass loss rate which 
is lower by a factor of ten. We investigated whether this 
could be due to the difference in the dust properties used. 
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Fig. 1. A fit to the spectral energy distribution of W Hya (solid line) with a three-dust component of amorphous 
silicates, AI2O3 and MgFeO. Open circles are from simbad (V magnitudes) and ground based observations (Wilson 
et al. I1972[l while filled squares are IRAS PSC and crosses are SWS and LWS data. Note that the model does not 
account for the 13/mi emission feature. 



Zubko & Elitzur used the silicate optical constants derived 
for AGN dust (Laor & Draine lT993[) . which has its infrared 
characteristics based upon astronomical silicates (Drainc 
& Lee ll984|) . We have also analysed the observed spectra 
using dust parameters for the astronomical silicates and 
for Mgo.sFeo.sSiOs from the University of Jena database. 
The derived dust mass loss rates which fit the observed 
strength of the 10/xm dust feature agree with our value. 
Hence the difference in the derived dust mass loss rate is 
not due to the differences in the optical properties used. 
Rather, Zubko & Elitzur (2000) did not attempt to fit the 



observed 10/xm feature in detail, but using only ground 
based-photometry and the ISO-LWS spectrum. We em- 
phasize that the high mass loss rate of Zubko & Elitzur 
( 200CJJ) gives too strong 10 and 20/im silicate emission as 
shown in Fig. [21 We find that extending the calculation 
to include the dust envelope size observed by IRAS at 60 
and 100/im (Hawkins 1990J does not affect the far-IR flux 
as the density is far too low in the outer part to affect 
the SED. For our case, we fitted the silicate 10/im feature 
and the other dust species (except for the 13/^m emission 
feature) which givess a more reliable mass loss estimate 
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Fig. 2. The plot of a close-up of the dust features. The 
long dashed line represent a model with Zubko & Elitzur 
(2000) parameters. The heavy solid line is the three dust 
model (Fig. while the heavy dashed line is a model 
which includes spinel. The contribution from individual 
dust component is shown as dotted lines. From top to bot- 
tom : spinel, amorphous silicate, AI2O3, and Mgo.1Feo.9O. 

as dust mass loss rate is proportional to the strength of 
the 10/im silicate feature. However, since we did not fit the 
13/im feature, our estimate of the mass loss rate should be 
viewed as a lower limit, pending the correct identification 
of the feature. 

If we assume that the elemental abundance of the out- 
flow is that of the Sun, our estimated dust mass loss rate 
for each dust component indicates that aluminum con- 
denses out at a higher rate than the silicate, magnesium 
and iron : the latter three having abundances an order 
of magnitude higher than that of Al (Anders & Grevesse 
1989p . The amount of the silicate dust is the same as the 
aluminum dust while that of MgFeO dust is five times 
less (Table 2J. Thermodynamic consideration of the dust 
condensation sequence starts with gas phase condensation 
of AI2O3, leading finally to Mg- and Fe-incorporation into 
silicate grains fTielens ll990p . Alternatively, Al and Si con- 
dense out at the same time but as Al has a high affinity 
for oxygen, aluminum oxide will be formed at a higher 
rate than silicate (Stencel et al. I1990f) thereby giving rise 
to an 11/im feature which is as strong as the silicate fea- 
ture even though the abundance of Al is much lower than 
that of Si. In either case, the high fraction of Al in dust 
as compared to Si, reflects a freeze out of the condensa- 
tion process well before all the silicon has been consumed 
(e.g., well before magnesium-iron silicates condense out). 
W Hya is prototypical for a set of O-rich stars which have 
10-20/im spectra which are not dominated by the well- 



known broad (magnesium/iron) silicate emission features 
but rather by a set of oxide features fCami l2"UU2"|l . We sur- 
mise that this difference reflects a steep density gradient 
in the dust condensation layer in these sources as com- 
pared to typical Mira variables. The presence of such a 
steep density gradient is also indicated by the asymmetric 
light curve of these types of stars. 

3.2. The molecular envelope 

Using the total dust mass loss rate obtained by fitting 
the SED which drives the gas to the terminal velocity, we 
obtain a dust-to-gas mass ratio of 8.7 10~ 3 , giving a gas 
mass loss rate of 3.5 10 -8 M Q yr^ 1 (Model 1) from solving 
the equation of motion, assuming that the gas and dust 
are fully momentum coupled. 

(i) 



dv 



where v is the gas velocity, G is the gravitational con- 
stant, Af* is the mass of the star, and T is the ratio of 
the radiation pressure force on dust grains to the gravita- 
tional force. Following Justtanont et al. I|1994[) . T can be 
written to include the dust size distribution, n d (a), and 
Q{a, A), the radiation pressure efficiency for a grain size, 
a, at wavelength A as 



a d Q{a, A)L(A) n d (a) 



A-kcGM* 



dX da 



P 



(2) 



where a d is the geometric cross section of the dust, L(X) 
is the luminosity at the wavelength A, p is the gas density 
(see Justtanont et al. (1994) for full details). We use ra- 
diation pressure for purely silicate dust in the calculation 
of the equation of motion. As a matter of fact, we find 
that the dust velocity is much greater than that of the gas 
due to the drift velocity in stars with such a low mass loss 
rate. 



"drift (a) = 



Q(a,X)L(\)v 
Mc 



da 



(3) 



The drift velocity is calculated for each grain size and 
then weighted by the size distribution of a -3 5 . The ratio 
for the drift velocity for the largest and smallest grains 
(0.25 - 0.005/im) is 8.1. In our calculation, the heating is 
dominated by the largest grains (heating oc v 3 ,). This re- 
sults in a high heating rate, therefore the calculated gas 
kinetic temperature is very high. In the fitting of the SED, 
we derived an M/t^ust in which we can ignore the drift 
velocity. However, if we take the drift velocity into account 
when calculating the dust mass loss rate and then solve 
for both the equations of motion and energy balance in a 
self consistent way, we derive a mass loss rate of 1.5 10 -7 
M yr _1 (Model 2) and the calculated CO line fluxes are 
higher than those observed for the low J transitions (see 
Table |3J) . From running a number of models with vari- 
ous different densities and temperatures, we find that the 
line flux ratio of the CO 2-1/1-0 is a sensitive measure of 
the assumed mass loss rate while the temperature plays a 
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Table 1. Input parameters for modelling the dust and gas 
mass loss rates in W Hya 



Stellar lTlcLSS 


1 Mn 




Stellar effective temperature 


2500 K 




kj LC11CL1 IdUlLlo 


4 10 13 cm 




Outflow velocity 


6 km s _1 




Distance 


115 pc 




Dust condensation radius 


2 10 14 cm 




Dust mass loss rates 


1.5 10~ 10 M yr _1 


- silicates 




1.3 10" 10 M yr _1 


- A1 2 3 




2.5 10" 11 M yr _1 


- MgFeO 


Dust-to-gas mass ratio 


8.7 10" 3 




Gas mass loss rate 


3.5 10" 8 M Q yr" 1 





Table 2. A comparison of between the observed lines and 
those from the models using parameters by Barlow et al. 
(1996) and those derived in this paper. We list the avail- 
able antenna temperatures (Tmb), line intensities (I) and 
line fluxes (F) from various sources. 



Transition 


Observed 


Barlow 


Model 1 


Model 2 


Model 3 






et al. 








Tmb (1-0) 


0.3K 


4.7 


0.06 


0.4 


0.1 


I (1-0) 


1.5" 


55.0 


0.8 


4.6 


1.57 


T MB (2-1) 


2.3K 


20.9 


0.7 


5.7 


1.9 


I (2-1) 


22. b 


229.6 


7.5 


51.4 


19.3 


F(16-15) 


3.3E-20 C 


8.0E-20 


1.6E-20 


5.9E-20 


2.9E-20 


F(17-16) 


2.3E-20 C 


8.6E-20 


1.7E-20 


6.2E-20 


3.1E-20 



minor role For the observed ratio of 14.7, an upper limit 
for the mass loss rate is 8 10 -8 M© yr^ 1 (Model 3, with 
the temperature derived from solving the energy balance 
equation). This model gives the best description of the ob- 
served line fluxes from low to higher transitions. As a test, 
we calculated CO line fluxes using this mass loss rate and 
a simple temperature law T oc (r/ri)~ a where a ranges 
from 0.5 to 0.8. The resulting line fluxes are insensitive 
to the change of temperature at this low mass loss rate 
(< 10%). This mass loss rate is within a factor of two of 
our detailed calculation and puts a stringent upper limit 
on the mass loss rate for W Hya. We also calculated the 
expected CO line fluxes using parameters from Zubko & 
Elitzur (2000), including their temperature structure. The 
calculated line fluxes are higher than those by Barlow et 
al. I[1996[) . reflecting the higher mass loss rate. 

We calculate the peak main beam antenna tempera- 
tures for Model 1 to be 0.1K for CO J=l-0 observed with 
the IRAM telescope (cf. 0.3K from Bujarrabal et al. I1989|> 
and 0.7K for J=2-l observed with a clear double-peaked 
profile (cf. a flat-top profile with a peak temperature of 
2.3K from Cernicharo et al. I|1997[l ). with a canonical value 
of the CO abundance of 3 10~ 4 . Our calculated gas kinetic 
temperature, obtained from solving the energy balance, in 
the outer part of the shell is very high due to the large drift 
velocity, hence the populations of low-lying levels are ex- 
cited to higher levels. This, coupled with a low density, 
results in very low antenna temperatures for J=l-0 and 2- 
1 lines. In order to reconcile these with the observations, 
we need to increase the mass loss rate to 8 10 -8 M© yr _1 
(Model 3). We also calculated expected line fluxes using 
input parameters from Barlow et al. I|1996|l and their tem- 
perature, and compared them to the IRAM results plus 
J=16-15 and J=17-16 in their ISO-LWS spectrum. The 
results of these calculations are summarized in Table [2] 
Our model predicts line fluxes which are lower than the 
observed values, which is due to the low derived mass loss 
rate for lower transitions while Barlow et al.'s parameters 
predict too high line fluxes by a factor of ten. However, we 
find that the high rotational lines are insensitive to the in- 
put mass loss rate. The fact that there is a combination of 
different types of dust, not just silicates, may lead to dif- 



a CO line intensity (K km s *) from Bujarrabal et al. H1989H 
b CO line intensity (K km s" 1 ) Cernicharo et al. (119971 
c ISO : Barlow et al. itT^oTl flux in W cm" 2 

ferent radiation pressure. A change in this parameter will 
affect the drift velocity and hence the heating rate and the 
gas kinetic temperature. Another factor which affects the 
CO line fluxes is the adopted turbulent velocity (Kemper 
et al. l2003f) . We assume a constant value of 1 km s _1 in our 
calculation but in reality this could vary within the cir- 
cumstellar envelope. An important issue which we touched 
upon in the previous section is the dust drift velocity. We 
calculated the drift velocity of ~ 15 km s _1 for the av- 
eraged radiation pressure for the grains (< Q(a, A) > = 
0.0057). However, the largest grains with a high drift ve- 
locity may no longer be momentum-coupled to the gas 
(fdrift ~ 100 km s~M). Hence these grains will not partake 
in the heating of the gas, and hence the temperature struc- 
ture of the envelope should be steeper. Further progress 
in this area has to await a more detailed radiative transfer 
calculation that includes a realistic model for the heating 
and cooling. 

Our value of the gas mass loss rate agrees within a 
factor of two to three with those from Wannier & Sahai 
(1986J and Olofsson et al. (j2002fl after scaling for the 
distance and velocity assumed here. However, our value 
falls short by an order of magnitude compared with that 
from Barlow et al. (1996) who calculated radiative trans- 
fer for H2O lines in the LWS wavelength range and from 
Neufeld et al. ill 996; who modelled H2O lines from the 
SWS. The gas mass loss rate derived recently by modelling 
H 2 emission lines of W Hya by Zubko & Elitzur QXfflty 
is slightly higher than that from Barlow et al. (1996). In 
an attempt to reconcile these differences, we can say that 
our assumed CO abundance of 3 10 -4 may be too high for 
this object. But we need to reduce this by a factor of 30 
to bring the derived mass loss rate in line with those from 
the water lines. However, the column density of the CO 
seen in absorption in SWS seems to contradict the very 
low abundance (Table |3J). The abundance of water used 
in these calculations may be too low. By increasing the 
water abundance and lowering the mass loss rate in those 
cases, or by increasing our mass loss rate and lowering 
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the CO abundance in our case, a convergence in mass loss 
rate can be found between these approaches. Also, it is 
difficult to get the kinematic information from ISO LWS 
observations of H 2 lines because they are not resolved. 
Although SWAS and Odin observations of the 557GHz 
line should be resolved, currently the signal-to-noise ratio 
is not good enough to be confident about the line pro- 
file. In the future it will be possible to study water lines 
in more detail using the heterodyne instrument (HIFI) 
aboard the Herschel Telescope due to be launched in 2007 
which has a resolution of 1 km s _1 or better. Equally im- 
portant, the collisional cross sections between ortho- and 
para-H^O and ortho- and para-H2 for high temperature 
extending to high quantum numbers will be available as 
part of the study project in conjunction with HIFI. In the 
case where the lines are subthermally excited, the input 
collisional rates have a strong influence in the calculation 
of radiative transfer. 



4. Modelling absorption bands 

The near-IR spectrum of W Hya shows several absorp- 
tion bands due to various molecules, such as CO, H2O. 
In order to see these bands clearly, we divided the spec- 
trum by our fit to the SED and normalised it at 4/xm. We 
have modelled these absorption bands assuming thermo- 
dynamic equilibrium. The molecular data used are taken 
from Lc Floch (1991) for the CO molecule, Schwenke 
(1998, private communication) for OH and H2O molecules 
and Hitran (Rothman et al., 1998) for others. The ab- 
sorption line profiles are convolved with a Voigt profile 
with the appropriate full width at half maximum for di- 
rect comparison with the observation (see Helmich[l996). 
From our SED fit to the spectrum, we can define the stel- 
lar continuum satisfactorily up to A < 9/xm. By fitting 
the observed continuum-divided spectrum with our mod- 
els, we can derive an excitation temperature (T x ) for each 
molecule, along with the column density, N, assuming a 
Doppler parameter, bp. In this paper, we try fitting the 
different molecular bands independently and derive the 
best set of parameters for each band. A more sophisti- 
cated approach has been made by Cami l|20()2[) using a 
chi-square technique. 

We fitted, in all, six different molecules from 2.3-9/im. 
The input parameters are listed in Table while the fits 
are shown in Figs. |21 and From the temperatures, we 
can deduce that these molecules originate from different 
layers in the circumstellar shell. We see bandheads of CO 
Av=2 from u=5-3 up to v=8-6 which indicates a very 
hot gas. This is reflected in the input parameter for CO 
which has a temperature of 3 000K. The column density 
for the band is high, 10 21 cm -2 also indicating that CO 
bands originate from a hot, dense region very close to the 
stellar photosphere. We believe that the gas is possibly 
excited by shocks since the effective temperature of the 
star is below the excitation temperature which gives rise 
to the observed v=8-6 band. This is also consistent with 
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2.5 3 3.5 4 

Wavelength ((jm) 

Fig. 3. ISO SWS observation in band 1 (dotted line) 
compared with the model absorption of CO Au=2, H2O 
stretching modes and OH absorption bands (solid line). 
Each absorption component is also depicted separately. 



the inferred Doppler width of 9 km s 1 used to model the 
CO absorption band. 



1 111 




Wavelength (yU.m) 



Fig. 4. ISO SWS observation in band 2 (dotted line) com- 
pared with the model absorption of C0 2 , H 2 bending 
mode, SO2 stretching mode and SiO v—1-0 band (solid 
line). Each absorption component is individually shown. 
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Table 3. Parameters for the spectral absorption calcula- 
tion. For each molecular band, we also indicate the wave- 
length coverage. 





\ 

A 


T 

-t X 


N 


VD 




(/im) 




\ cra ) 


(km s _1 ) 


CO Av=2 


2.3-2.6 


3 000 


1 10 21 


9 


H2O Vl + IAj 


2.3-4.0 


1000 


510 20 


7 


OH 


2.7-4.0 


1500 


710 17 


8 


C0 2 V3 


4.0-4.6 


1000 


810 17 


7 


CO v=l-0 


4.2-5.5 


2 000 


510 18 


8 


H 2 V 2 


4.0-9.0 


300 


510 18 


3 


SO2 


7.1-7.7 


500 


210 17 


3 


SiO u=l-0 


7.2-9.0 


1000 


no 19 


7 


C0 2 u 2 


14.8-15.0 


1000 


no 17 


7 



We also see the stretching modes of H 2 molecules 
at 3/xm. We modelled these bands assuming an ortho to 
para ratio of 3. We are able to reproduce the absorption 
band relatively well with a gas temperature of 1 000K for 
the H2O stretching modes. However, the band is narrower 
than the observed one, especially in the red wing. The 
water band overlaps the OH ro- vibrational bands ^3/2 1 — 
n 3 / 2 and — n 1 / 2 0) which cover from 2.9 to 4//m. 

The highly excited lines of OH also indicate that they 
come from a warm region, close to the star where H 2 
is probably photodissociated by shocks. We also modelled 
the C0 2 stretching mode at 4.25/im. The band seems to 
come from the same region as the others. The width of 
the band is consistent with a gas of 1 000K. 

The fit to the bending mode of H2O at 6/im is, how- 
ever, not consistent with the observed data. In view of 
the various ratios of adjacent lines, we estimate that the 
excitation temperature is only 300K for a column density 
of 510 18 cm -2 . The observed absorption is much deeper 
than calculated in the blue wing where the CO fundamen- 
tal band is present. The 4.6/im CO R-branch blends with 
the CO2 absorption wing while its P-branch blends with 
H2O. At the resolution of the observations, individual ab- 
sorption lines are not resolved. We roughly estimate the 
excitation temperature of the CO fundamental absorption 
to be ~ 2 000K and put an upper limit for the column den- 
sity of 10 19 cm -2 . From very high resolution studies of CO 
lines in Mira variables (Hinkle et al. 119821) . different bands 
come from different regions, with the fundamental band 
coming from a cool region while the overtone bands come 
from a hotter region, close to or from the pulsating photo- 
sphere. We can also fit the SiO band at 8/im with a warm 
gas of 1 000K, possibly from the same layer as the H2O 
stretching modes. This is consistent with emission coming 
from gas interior to the dust forming region where SiO is 
not bound in dust grains. There is evidence for the weak 
SO2 absorption band at 7.3/im. This band was discovered 
in a few O-rich AGB stars by Yamamura et al. (|1999ll . We 
also see CO2 ^2 absorption band at 15/^m which we will 
discuss in more detail in the next section. 



Table 4. Line fluxes of rotational water lines detected in 
the ISO-SWS wavelength range. 



wavelength 
( //ml 


line flux 
(W cm -2 ) 


transition 




32.98 


1.9E-18 




^bl 


35.45 


5.3E-19 


^33 ~U4 , 1 43 


— fvj/i 

^34 


35.92 


4.3E-19 


652 — 541, 651 


— 542 


39.35 


1.6E-18 


5so — 44i , 551 


— 440 , 642 ~ 533 


40.31 


4.8E-19 


643 — 532 




40.71 


1.4E-18 


432 — 3o3, 633 


— 524 


43.91 


6.6E-19 


541 — 4 3 2 




44.14 


6.1E-19 


542 — 431 




44.68 


4.4E-19 


836 — 725 




45.10 


6.8E-19 


523 - 4i4 





Although the overall fit to each band is good, there is 
still some residue which we are not currently able to ac- 
count for. There is an absorption band just longward of the 
4.25/im CO2 band. It is likely that the hotter component of 
water from the stretching modes contributes between 4.5- 
6/im. We compared our results to that from Cami (j5002 
Table 6.2) and found that although our derived excitation 
temperatures agree within a factor of two, the column den- 
sities are very different (up to two magnitudes difference) 
in the cases of OH and SiO, while others agree reason- 
ably well. This may lie in the fact that we adopt different 
Doppler broadening for different regions, i.e., the hotter 
region has a higher thermal broadening than the cooler 
one, while Cami (2002) assumed a constant velocity of 3 
km s _1 throughout. As it turns out, different combinations 
of column density and Doppler width can give similar re- 
sulting absorption. The depth of the absorption is very 
sensitive to the adopted Doppler velocity used, more so 
than the column density. 

5. Molecular emission lines 

Both Barlow et al. <fiT?9l)j) and Neufeld et al. ifEMjl re- 
ported emission due to H2O rotational lines from LWS 
and SWS FP mode observations, respectively. From the 
full SWS scans, we detected 10 more H2O rotational lines 
in the region between 30-45/im. We estimated the line 
fluxes by assuming a gaussian profile and these are re- 
ported in Table 01 Some lines are a blend of the ortho- 
and para-water as the wavelengths of their transitions are 
too close to be resolved by the SWS scan. Barlow et al. 
(1996) also noted some CO lines present in their LWS 
spectrum. Surprisingly, only two CO lines have been de- 
tected in LWS (J=16-15 and J=17-16, see Sect. 3). In this 
section, we will limit our discussion mainly to new results 
on C0 2 . 

5.1. The SWS grating spectrum of CO2 

We detected CO2 in W Hya between 13-16/J,m from the 
full scan SWS AOT1. We see these C0 2 bands in all of 
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14.98 



Fig. 6. SWS-FP observations of CO2 in W Hya with the wavelength corrected for the heliocentric velocity. The 
top panels show the partly resolved Q-branches of the 13.87 and 16.18/im bands and the bottom panel shows the 
14.97/xm. band. Flux units are 10~ n W m~ 2 /im" 1 . Vertical lines indicate positions of Q lines. The first few Q 
lines are not resolved. 



the stars which exhibit the so called 13/im dust feature in 
our guaranteed time program (PI. de Jong) . From this we 
proposed that objects with the 13/^m dust have CO2 asso- 
ciated with it (see Justtanont et al. 1998, Fig. 1). On this 
basis, we obtained further observing time to scan objects 
with this dust feature at the highest grating resolution 
(AOT6). All the objects observed in the supplementary 
program confirm our idea, i.e., all stars with the 13/im 
dust show CO2 bands. However, the reverse is not always 
the case, e.g., one of the stars in our program, o Ceti, has 



no 13/im dust but shows CO2 emission and IRC+ 10011 
also shows CO2 while the 10/im silicate is self-absorbed 
(Markwick & Millar I207)0"| . Cami et al. pW)| show a de- 
tailed calculation of CO2 emission bands in EP Aqr, one 
of the stars in our program. If the association between 
the 13/zm dust feature and CO2 gas is real then CO2 is 
very prevalent because Sloan et al. i|199ri|l argue that up to 
50% of O-rich semiregular variables show such a dust fea- 
ture. The formation of CO2 is thought to be via CO+OH 
(Nercessian et al. 115551 Willacy & Millar IT59T|) where OH 
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12 13 14 15 16 



Wavelength (/im) 

Fig. 5. The full grating resolution spectrum of W Hya, 
showing the 13/^m dust and CO2 bands. The insert show 
the modelled absorption (solid line) compared to the ob- 
served spectrum (dotted line). 
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Fig. 7. The rotation diagram of CO2. Crosses are from the 
01 1 0-00°0, triangles represent the 02°0-01 1 and squares 
represent 10°0-01 1 0. The solid lines are the best fit for 
these bands. 



is the product of either photodissociation or of shock de- 
struction of H2O molecules. 

In W Hya, the transitions 10°0-01 1 at 13.87/xm and 
02°0-01 1 at 16.18/mi are in emission while 01 1 0-00°0 at 
14.97/im is in absorption. We modelled the absorption 
band which is shifted from the normal position as seen 
in emission in other stars. We attribute this shift towards 
a shorter wavelength to the importance of highly excited 
bands. To excite these bands (e.g., 02 2 0-01 1 0, 03 3 0-02 2 0, 
04 4 0-03 3 0), we need the gas temperature to be ~ 1 000K 
(Fig. [SJ). These bands have their Q-branches at succes- 
sively shorter wavelengths than the 01 1 0-00°0 band. By 
including these bands in our calculation, we are able to 
match the peak absorption at 14.957/im as observed. One 
thing to note from our observations of objects with CO2 
is that only stars with the 14.97/im absorption show the 
corresponding 4.25/im absorption. Those which show the 
14.97/im in emission lack the stretching mode counterpart 
(for the full sample, see Cami l|2l)U2[) ). This probably in- 
dicates that the radiative pumping by the 4.25/im band 
plays a very minor role, if at all, in populating the levels 
which result in the 13-16/im emission. 

5.2. The Fabry-Perot observation of CO2 bands 

We subsequently obtained high resolution spectra of these 
CO2 bands using the SWS-FP which resolved the individ- 
ual Q branches (Fig.|HJ). The small dip at 13.877/im in the 
band 10°0-01 1 in this figure is an instrumental effect due 
to the end of the scan of FP wavelength setting. The re- 



solved Q-lines are listed in TableEJwhere the wavelength, 
the energy of the upper level, the measured line flux and 
the identification are given. The flux errors given repre- 
sent the estimated uncertainties in the continuum and do 
not take into account possible calibration errors so that 
they should be taken as lower limits. 

It is quite remarkable though that we can discern emis- 
sion of the Q-branch lines in the depth of the 14.97/im 
band while the absorption blends into a single band. As 
discussed in Sect. 5.1, this absorption is formed in a hot 
(1 000K) gas layer in the outer photosphere and lines may 
be smeared out by turbulence (6d = 7 km s _1 , see TableOJ 
and/or optical depth effects. The first of the individual On- 
line to be resolved in emission is Q10. The close spacing 
between Q2-Q8 means that even at this resolution we can- 
not separate them from each other. We conclude that the 
emission lines are produced in a gas layer which is more 
extended and lying farther out in the circumstellar shell 
than the hot layer which is responsible for the absorption 
band at 14. 95-14. 98/im, in line with earlier suggestions by 
Justtanont et al. I|1998|) and Cami et al. H2000[> . 

We have Fabry-Perot observations of individual Q- 
branch lines for three bands of CO2 (see FigEJ. In the 
fundamental band we detected lines up to Q42 and in the 
other bands up to Q24. We can use the observed flux of 
the individual line in each band to construct a rotation 
diagram from which we can derive a rotational excitation 
temperature T rot characterizing the population distribu- 
tion within the upper vibrational level of each band, using 
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Table 5. A list of the CO2 emission lines detected from 
the FP observations. 



1 //mi 


E 

(K) 


linp fliTx - 

f10~ 19 W rm" 2 ) 


transition 




10.0 1 JU^ 


9059 0606 


2.6±0.5 


1 n u n 

1U 


ni A n 

Ul u 


01 n 

IU 


1 3 87658 


9084 8831 


1.0±0.5 


1 n°n 

1U 


ni 1 o 

Ul u 


019 

1^ 


1 R77C\R 

-LO.O 1 ( DO 


911^1 0^4 


2.2±0.2 


1 n°n 


ni 1 n 

Ul u 


D14 


1 S R7RQA 
10.0 1 oy^ 


91 49 QQ74 


2.7±0.2 


1 n°n 


m x n 

Ul u 


Ol 6 

IU 


IO.OOUOU 


91 &Q 9889 
zioy.ZiOo^i 


1.9±0.2 


1 n°n 

iu u- 


m x n 

Ul u 


ni s 


1 3 RR~\ Q4 


OOQQ f)fi77 


^ 3+1 Q 


1 n°n 

iu u- 


m 1 n 

Ul u 


n9n 


1 ^ RR1C\7 

IO.OOOU ( 


OOQl C= 1 
ZtZtO 1- .OOO -L 


3 1 -M 5 


1 n°n 

iu 


ni 1 o 

Ul u 


D99 


l-O .oooov 


Zi / -± . yj y u 


4.3±1.1 


1 n°n 

iu u- 


ni 1 n 

Ul u 


D94 


10.00 ( uu 


Zioyi.oouj 


1 R+n 5 

i .uznu.u 


1 n°n 


ni 1 n 

Ul u 




1 4 Q499Q 


1 Q7fi 941 ^ 


1 fi+n 3 

i .uznu.o 


m 1 n 


nn°n 

uu u 


049 


1 4 Q4R0Q 


1 889 Q0fS4 
ioo//.yuu4 




m x n 

ui u- 


nn°n 

-uu u 


04n 


1 4 Q4Q79 


1 7Q4 fl^Q4 


2.1±0.3 


ni 1 n 

ui u- 


nn°n 

-uu u 




1 4 Q^l 7 


1 70Q 701 R 


1 c_i_n q 
_l .vjznu.o 


ni 1 n 


nn°n 

uu u 






1 fi9Q R^A 


n Q+D 4 


m 1 n 

ui u- 


nn"n 

uu u 


034 


14:. £70 iJ'J^J 


1 ^4 4fi1 1 


1.0±0.3 


m x n 


nn°n 

uu u 


039 


1 4 Qfi94^ 


1400.00UO 


1 Q+n ^ 


m x n 


nn°n 

uu u 


o^n 




1 41 7 1 Q4^ 


2.1±0.4 


ni 1 n 

ui u- 


nn°n 

-uu u 


098 


1 4 Q(S7(SQ 

1^1. C/U 1 \JU 


1 355 3043 


2.1±0.5 


ni 1 n 

Ul u 


nn°n 

uu u 


096 


1 4 Q700^ 


1 9Q7 Q1 0Q 


1.2±0.4 


m 1 n 

Ul u 


nn°n 

uu u 


D94 


1 4 Q799'} 


1 94^ 01 ^ 


1.2±0.4 


ni 1 n 

Ul u 


nn°n 

uu u 


D99 


1 4 97499 


1 1 96 61 87 


1.3±0.5 


ni x n 

Ul u 


nn°n 

uu u 


09n 


1 4 97603 


1 1 59 791 6 


2.1±0.7 


ni x n 

Ul u 


nn°n 

uu u 


Ol 8 
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1 4 0778^ 


1 1 1 Q ^94(S 

_1_ _L _1_ O. OZr4:U 


9 fi+fl 7 


ni x n 

Ul u 


nn°n 

uu u 


Ol fi 


1 4 Q7Q0Q 


1 078 4985 

\-\J 1 0. L ±£i<JO 


3.8±0.3 


O^O- 


-00°0 


014 


1 4 98034 


1 048 033Q 


2.2±0.5 


o^o- 


■00°0 


019 


14.98141 


1099 1414 


2.7±0.7 


Ol^- 


-00°0 


oin 


16.18366 


1937.0742 


4.2±0.2 


02°0- 


-Ol^ 


Q12 


16.18476 


1967.4081 


2.2±0.2 


02°0- 


•01*0 


Q14 


16.18603 


2002.2345 


1.7±0.2 


02°0- 


-Ol^ 


Q16 


16.18747 


2041.5522 


2.5±0.2 


02°0- 


•Ol^ 


Q18 


16.18908 


2085.3610 


6.0±2.0 


02°0- 


•Ol^ 


Q20 


16.19087 


2133.6598 


4.1±0.5 


02°0- 


•Ol^ 


Q22 


16.19283 


2186.4477 


5.6±0.1 


02°0- 


-Ol^ 


Q24 



the relation 

E 3 /kT mt = ln(N v /$(T Iot )d 2 ) - In^izFji/Ajihvjigj) 

= ln{N v /^{T mt )d 2 ) - Yj (4) 

where Ej and gj are the energy and statistical weight of 
the upper state, Fji and Vji are the line flux and frequency 
of the transition, Aji is the spontaneous transition prob- 
ability, N v is the total number of molecules in the upper 
vibrational level v, $(T rot ) is the partition function, and 
d is the distance to W Hya. Fig. shows the rotation di- 
agram of all the resolved Q-lines. There is considerable 
scatter in the data points due to uncertainties in the cal- 
ibration (e.g., joining of different pieces of the spectrum 
within each band) and in the continuum. The offset be- 
tween the data points in the 01 1 level and in the 02°0 and 
10°0 levels is due to differences in the number of molecules 
in these levels. 

In Table we collect a number of quantities derived 
from the fits in Fig.0 The derived rotation temperatures 
are 435 K for the 01 1 level and 460 K for the other two 



Table 6. Quantities derived from the rotation diagram 



Level 


T rot (K) 


$(T rot ) 


Y 


Nv 


10 u 


460 


4.73 


-7.57 


3.0 10 as 


02°0 


460 


6.56 


-7.50 


4.2 10 3S 


01*0 


435 


47.2 


-3.20 


2.4 10 41 



levels combined. The partition functions for each vibra- 
tional level are computed assuming that the levels are 
populated in thermal equilibrium at the rotational tem- 
perature. The quantities Yo are the intercepts of the fits 
in Fig. □ at Ej = 0. Putting Ej =0 and Y,- = Y in Eq. 
(3) we derive values of N„, the total number of molecules 
in each vibrational level, listed in column 5 of Table |5J 
These numbers are a lower limit as we assumed that the 
lines are optically thin. 

From these, we then find the vibrational excita- 
tion temperatures T vib (10°0-01 1 0) = 150 K and T vib 
(02°0-01 1 0) = 170 K, using the relation T vib = 1.44 
E y /ln(N 3 /N I ) withE(10°0-01 1 0) = 720 cm" 1 andE(02°0- 
01 1 0) = 617 cm -1 and taking into account that the num- 
ber of rotational levels in level 01 1 is twice as large (J 
= even and odd) as in levels 10°0 and 02°0 (J = even 
only). Assuming that the relative number of molecules in 
the ground vibrational state 00°0 and the first excited vi- 
brational state 01 1 is also characterized by T v ib ~ 160 K 
we find a total number of CO2 molecules of 5 1 43 . 

Clearly the temperature characterizing the distribu- 
tion of rotational levels in the vibrational states of CO2 
(~ 450 K) is different from the temperature characterizing 
the distribution of molecules over the vibrational states 
(~ 160 K) so that the CO2 molecules are close to, but 
not in thermal equilibrium. This non-LTE behaviour can 
be understood on the basis of the symmetry of the CO2 
molecule which does not allow pure rotational radiative 
transitions and the magnitudes of the cross sections for 
(de-)excitation of CO2 by collisions with H2, which are 
about three orders of magnitude larger for pure rotational 
transitions than for ro-vibrational transitions (Banks & 
Clary I*lW7|) . It can be shown by more detailed non-LTE 
multi-level line transfer calculations (de Jong, in prepara- 
tion) that at a few stellar radii from the star and at H2 
densities above ~ 10 8 cm -3 the rotational levels are more 
strongly coupled to the gas kinetic temperature (~ 400 K) 
while the population of the vibrational levels is dominated 
by radiation trapping in the ro-vibrational lines (the ones 
that we observe) resulting in a vibrational temperature 
somewhat below the kinetic temperature. 

We suggest that CO2 absorption and emission origi- 
nate in two different regions of an extended dynamical at- 
mosphere : one close to the star in the outer photosphere 
where the temperature is 1 000K, giving rise to the absorp- 
tion seen at 4.25 and 14.95/im, and the other farther out 
at few stellar radii, where the density is still high enough 
to thermally distribute the pure rotational levels at ~ 450 
K. Assuming that the cooler CO2 layer is located at ~ 3 
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stellar radii we find, from the total number of molecules 
observed a column density N(CC>2) ~ 10 15 two orders of 
magnitude smaller than that derived for the hot CO2 gas 
(see Table 3) and also several orders of magnitude smaller 
than found for the warm gas in the star EP Aqr by Cami 

et al. (T2UUUjl 

As proposed earlier by Tsuji et al. I|1997f) . Yamamura 
et al. l|1999|l . and Cami et al. (2000]), we believe that the 
emission line gas detected here is part of the warm gas 
layer located between the stellar photosphere and the dust 
condensation zone. This gas is probably related to the 
blown-up atmosphere observed by Reid & Menten (1990) 
and found by Homer et al. I|1998|l in their calculations of 
the dynamical atmospheres of pulsating stars. How C0 2 
molecules are formed in this gas is unclear but they may 
well be a left-over from the preceding photosphere phase. 
According to Hofner et al. I)1998[l . this gas is expected to 
have turbulent velocities of order 1 km s . Eventually a 
little further out (at ~ 5 stellar radii) the dust is formed 
which then accelerates the flow to 5-10 km s _1 . The CO2 
is probably quickly photodissociated in the expanding en- 
velope. We note here that the line widths of the emis- 
sion are much narrower than those of the absorption, sug- 
gesting that the emission comes from a cooler (since we 
do not see higher overtone bands in emission), much less 
turbulent layer than that which produces the absorption. 
Whether the hot and cool layers of CO2 are manifesta- 
tions of the same gas needs to be investigated further in 
a detailed model for CO2 excitation. 

6. Summary 

By fitting the SED of W Hya from the near- to far-IR, we 
derive a dust mass loss rate of 310 -10 M Q yr _1 using a 
combination of three different dust species of amorphous 
silicate, AI2O3 and MgFeO. We also tried adding a fourth 
component of spinel to account for the 13/xm dust fea- 
ture but in order to fit it, we also expect to see another 
strong peak at Wfim which is not present in the observed 
spectrum. From this, we conclude that spinel cannot be 
the carrier of the 13^m dust emission feature. Hence, the 
debate on the carrier of the 13fim feature is still open. 
With the estimated value of the our dust mass loss rate, 
we get a gas mass loss rate range of (3.5-8) 10~ 8 M Q yr _1 
from simple consideration of momentum dust driven wind. 
Our derived mass loss rate agrees well with those previ- 
ously derived to explain the CO emission, taking into ac- 
count the differences in the adopted distance and velocity. 
However, the discrepancy still exists between these mass 
loss rates and those derived from the H2O emission, the 
latter being consistently higher. The difference in derived 
mass loss rates from CO and H2O may simply be due to 
the abundances assumed for each species. 

We obtained the full scan SWS spectrum of W Hya in 
which we modelled the absorption bands due to various 
molecules. The derived excitation temperatures of differ- 
ent bands suggest that they originate from different layers 
of gas, from 300K to 3 000K. However, some molecules ex- 



ist over a wide range of temperatures and cannot really 
be stratified within a particular zone, particularly those 
which show multiple absorption bands. The exception may 
be OH which is a product of water molecules which are 
dissociated by shocks in the inner part of the envelope. 

We also observed CO2 bands in emission using SWS 
FP and were able to resolve the three strongest Q branches 
in the mid-IR. The resulting rotation diagram yields an ex- 
citation temperature of ~450K while the CO2 absorption 
at 14.97/im requires a gas of 1 000K to explain the shift 
of the peak absorption bluewards. The stretching mode at 
4.25/im also requires a layer of hot gas to explain its width. 
We find from simple calculation based on the rotation di- 
agram that the vibrationally excited levels may not be in 
LTE. Whether the CO2 absorption and emission lines can 
indeed be explained by temperature stratification in an 
extended dynamical atmosphere needs to be investigated 
in a more detailed radiation transfer model. 
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